Author: Ranjit Deshmukh 11/13/2016 This script does the following 1) Creates the results table using outputs from the results csv from the Renewable Energy Value model

## [1] "/Users/ranjitster/Dropbox/renewable_energy_value/renewable_energy_value"

Define input files…

Specify scenarios to process…

User parameter inputs…

Function to grab legend…

Import results and other parameter from csv files…

Creating empty results data tables and lists for plots…

Capital recovery factor…

Preparing the results input table…

Variables description

VRE stands for variable renewable energy.

Energy generation, capacity, and dispatch costs - inputs from REV model script

dispatch_cost: Total annual dispatch cost
ann_gen_total_MWh: Annual total energy generation (includes losses due to battery storage)
ann_gen_vre_nocurt_MWh: Uncurtailed VRE generation potential
ann_gen_vre_MWh: Annual VRE energy generation
ann_curt_vre_MWh: Annual curtailed VRE generation
ann_gen_solarPV_MWh: Annual Solar PV energy generation
ann_gen_wind_MWh: Annual Wind energy generation
ann_gen_coal_MWh: Annual Coal energy generation
ann_gen_gas_ccgt_MWh: Annual CCGT energy generation
ann_gen_gas_ct_MWh: Annual CT energy generation
ann_gen_diesel_MWh: Annual Diesel energy generation
ann_gen_hydro_MWh: Annual Hydro (including storage hydro and RoR) energy generation
ann_gen_nuclear_MWh: Annual Nuclear energy generation
ann_gen_other_MWh: Annual ‘Other’ energy generation including biomass, bagasse.
ann_discharge_bat_storage_MWh: Annual battery storage discharge
ann_charge_bat_storage_MWh: Annual battery storage charge
new_capacity_coal_MW: New coal capacity build
new_capacity_gas_ccgt_MW: New ccgt capacity build
new_capacity_gas_ct_MW: New ct capacity build
capacity_vre_MW: VRE capacity
capacity_solarPV_MW: Solar PV capacity
capacity_wind_MW: Wind capacity
capacity_bat_storage_MW: Battery storage capacity in MW
capacity_bat_storage_MWh: Battery storage capacity in MWh (MWh/MW gives number of hours of battery storage)

Total annual costs

ann_cap_cost_coal_mil: Annual capital cost of coal in millions ann_cap_cost_ccgt_mil: Annual capital cost of CCGT in millions ann_cap_cost_ct_mil: Annual capital cost of CT in millions ann_cap_cost_solar_mil: Annual capital cost of solar PV in millions ann_cap_cost_wind_mil: Annual capital cost of wind in millions ann_cap_cost_battery_mil: Annual capital cost of battery storage in millions total_cost_new_conv: Total annual cost of new conventional capacity (annualized cost)
total_cost_vre: Total annual cost of VRE in million USD (annualized cost)
total_cost_capacity: Total cost of conventional + VRE capacity i.e. total_cost_new_conv + total_cost_vre
total_cost_energy: Total annual cost of energy generation from all generation sources (result from economic dispatch model)
total_cost: Total annual cost i.e. total_cost_capacity + total_cost_energy (annualized fixed costs of existing conventional capacity is not included here. Cancels our for additional costs of VRE)
system_cost_vre: Total annual additional cost of vre i.e. (total cost of VRE scenario - total cost of zero VRE scenario)

Costs per curtailed or uncrutailed VRE

cost_vre_pMWh: Cost of VRE per MWh of VRE after curtailment i.e. (total_cost_vre/ann_gen_vre_MWh)
cost_vre_nocurt_pMWh: Cost of VRE per MWh of VRE before curtailment i.e. (total_cost_vre/ann_gen_vre_nocurt_MWh)
cost_battery_pMWh: Cost of battery storage per MWh of VRE after curtailment i.e. (ann_cap_cost_battery_mil10^6/ann_gen_vre_MWh)
cost_battery_nocurt_pMWh: Cost of battery storage per MWh of VRE before curtailment i.e. (ann_cap_cost_battery_mil
10^6/ann_gen_vre_nocurt_MWh)
cost_vre_battery_pMWh: Cost of VRE and battery storage per MWh of VRE after curtailment i.e. (cost_vre_pMWh + cost_battery_pMWh) <– USE THIS FOR EITHER VRE OR BOTH VRE AND BATTERY
cost_vre_battery_nocurt_pMWh: Cost of VRE and battery storage per MWh of VRE before curtailment i.e. (cost_vre_nocurt_pMWh + cost_battery_nocurt_pMWh) <– USE THIS FOR EITHER VRE OR BOTH VRE AND BATTERY

Capacity and energy value

capacity_value_pMWh: Capacity value of VRE target per VRE generation after curtailment i.e. difference between annualized costs of new conventional capacity of VRE scenario and zero VRE scenario per MWh of VRE generation
energy_value_pMWh: Energy value of VRE target per VRE generation after curtailment i.e. difference between total annual energy costs of VRE scenario and zero VRE scenario per MWh of VRE generation
total_value_pMWh: Total of capacity and energy value of VRE i.e. (capacity_value_pMWh + energy_value_pMWh)

Additional system cost due to VRE

system_cost_vre_pMWh: Additional system cost due to VRE for a VRE scenario compared to a zero VRE scenario per MWH of VRE i.e. (cost_vre_pMWh + cost_battery_pMWh - capacity_value_pMWh - energy_value_pMWh)
system_cost_add_pMWh_load: Additional system cost due to VRE for a VRE scenario compared to a zero VRE scenario per MWH of Load served. (Load served includes battery storage losses in battery scenarios. Additional system cost is divided by ann_gen_total_MWh.)
system_cost_add_perc: Additional overall cost of VRE as a % of no VRE system cost

Cost or benefit of mitigating CO2 emissions

total_emissions_co2_milTonnes: Total CO2 emissions in million tonnes
grid_emissions_factor_kg_kWh: Grid emissions factor
cost_emissions_reduction_pTonneCO2_avg: Average cost of CO2 reduction in USD per tonne of CO2
cost_emissions_reduction_pTonneCO2_marginal: Marginal cost of CO2 reduction in USD per tonne of CO2 (for every 200 GW target)

Processing results…

[1] “base”

[1] “coal_55mingen”

[1] “coal_ccgt_0mingen”

[1] “wind10LC”

[1] “wind20LC”

[1] “wind30LC”

[1] “solar10LC”

[1] “solar20LC”

[1] “solar30LC”

[1] “wind30LC_solar30LC”

[1] “wind120HH_solar1A”

[1] “battery15”

[1] “battery30”

[1] “battery60”

[1] “battery15B25LC”

[1] “battery30B25LC”

[1] “battery15B50LC”

[1] “battery30B50LC”

[1] “battery60B50LC”

[1] “high_cost_coal”

[1] “wind30LC_solar30LC_coalHC”

[1] “load_modD0M0_energyOnly”

[1] “load_modD50M0_energyOnly”

[1] “load_modD25M25_energyOnly”

Processed Results Table…

   scenario scenario_main dispatch_cost ann_gen_total_MWh
1:     S0W0          base  4.607645e+12        2424553182
2:   S0W200          base  3.532123e+12        2424553182
3:  S50W150          base  3.583863e+12        2424553182
4: S100W100          base  3.669521e+12        2424553182
5:  S150W50          base  3.817777e+12        2424553182
6:   S200W0          base  4.037051e+12        2424553182
   ann_gen_vre_nocurt_MWh ann_gen_vre_MWh ann_curt_vre_MWh
1:                      0               0                0
2:              477166828       477093308            73516
3:              454821772       454776174            45593
4:              427217839       426206260          1011581
5:              392945562       375118489         17827089
6:              348007539       278446347         69561203
   ann_gen_solarPV_MWh ann_gen_wind_MWh ann_gen_coal_MWh
1:                   0                0       2226209347
2:                   0        477093308       1763116977
3:            87798653        366977521       1783235771
4:           174521298        251684961       1807136233
5:           246842981        128275507       1849931247
6:           278446347                0       1945090909
   ann_gen_gas_ccgt_MWh ann_gen_gas_ct_MWh ann_gen_diesel_MWh
1:             37431382            1890048           187717.2
2:             23986378            1291599           229093.1
3:             25922219            1358091           426239.2
4:             30079658            1711104           585240.6
5:             37249322            2553021           866416.0
6:             38237836            2808021          1135382.3
   ann_gen_other_MWh ann_gen_hydro_MWh ann_gen_nuclear_MWh
1:              0.00         128787887            30046800
2:           1139.25         128787887            30046800
3:              0.00         128787887            30046800
4:              0.00         128787887            30046800
5:              0.00         128787887            30046800
6:              0.00         128787887            30046800
   ann_discharge_bat_storage_MWh ann_charge_bat_storage_MWh
1:                             0                          0
2:                             0                          0
3:                             0                          0
4:                             0                          0
5:                             0                          0
6:                             0                          0
   new_capacity_coal_MW new_capacity_gas_ccgt_MW new_capacity_gas_ct_MW
1:               151140                        0                  13400
2:               126720                        0                  28800
3:               118800                        0                  25000
4:               119460                        0                  27400
5:               126720                        0                  24200
6:               144540                        0                  11000
   capacity_vre_MW capacity_solarPV_MW capacity_wind_MW
1:               0                   0                0
2:          200005                   0           200005
3:          200025               50015           150010
4:          200102              100079           100023
5:          200039              150010            50029
6:          200004              200004                0
   capacity_bat_storage_MW energy_bat_storage_MWh scenario_build
1:                       0                      0           0 GW
2:                       0                      0         200 GW
3:                       0                      0         200 GW
4:                       0                      0         200 GW
5:                       0                      0         200 GW
6:                       0                      0         200 GW
   scenario_split rank           date_time scenario_econ_dispatch
1:            0-0    1 2018-01-17 00:22:52                ClcC70m
2:          0-100    2 2018-01-17 14:02:44                ClcC70m
3:          25-75    3 2018-01-17 16:51:26                ClcC70m
4:          50-50    4 2018-01-17 18:07:51                ClcC70m
5:          75-25    5 2018-01-17 19:30:53                ClcC70m
6:          100-0    6 2018-01-17 20:59:34                ClcC70m
   scenario_sensitivity ann_cap_cost_coal_mil ann_cap_cost_ccgt_mil
1:              ClcC70m              19006.02                     0
2:              ClcC70m              15935.17                     0
3:              ClcC70m              14939.23                     0
4:              ClcC70m              15022.22                     0
5:              ClcC70m              15935.17                     0
6:              ClcC70m              18176.06                     0
   ann_cap_cost_ct_mil ann_cap_cost_solar_mil ann_cap_cost_wind_mil
1:            873.4057                  0.000                 0.000
2:           1877.1705                  0.000             21878.861
3:           1629.4883               3933.600             16409.829
4:           1785.9191               7871.055             10941.673
5:           1577.3446              11798.049              5472.751
6:            716.9748              15729.997                 0.000
   ann_cap_cost_battery_mil total_cost_new_conv total_cost_vre
1:                        0            19879.42           0.00
2:                        0            17812.35       21878.86
3:                        0            16568.71       20343.43
4:                        0            16808.14       18812.73
5:                        0            17512.52       17270.80
6:                        0            18893.03       15730.00
   total_cost_capacity total_cost_energy total_cost system_cost_vre
1:            19879.42          70886.84   90766.26              NA
2:            39691.21          54340.36   94031.56        3265.301
3:            36912.14          55136.35   92048.49        1282.230
4:            35620.87          56454.17   92075.04        1308.773
5:            34783.32          58735.03   93518.35        2752.086
6:            34623.03          62108.48   96731.51        5965.246
   cost_vre_pMWh cost_vre_nocurt_pMWh cost_battery_pMWh
1:           NaN                  NaN               NaN
2:      45.85866             45.85160                 0
3:      44.73284             44.72835                 0
4:      44.13996             44.03544                 0
5:      46.04092             43.95214                 0
6:      56.49202             45.20016                 0
   cost_battery_nocurt_pMWh cost_vre_battery_pMWh
1:                      NaN                   NaN
2:                        0              45.85866
3:                        0              44.73284
4:                        0              44.13996
5:                        0              46.04092
6:                        0              56.49202
   cost_vre_battery_nocurt_pMWh capacity_value_pMWh energy_value_pMWh
1:                          NaN                  NA                NA
2:                     45.85160            4.332646          34.68186
3:                     44.72835            7.279860          34.63350
4:                     44.03544            7.206089          33.86312
5:                     43.95214            6.309745          32.39460
6:                     45.20016            3.542470          31.52623
   total_value_pMWh system_cost_vre_pMWh system_cost_add_pMWh_load
1:               NA                   NA                        NA
2:         39.01451             6.844155                 1.3467638
3:         41.91336             2.819476                 0.5288523
4:         41.06921             3.070751                 0.5397998
5:         38.70434             7.336577                 1.1350898
6:         35.06870            21.423322                 2.4603485
   system_cost_add_perc vre_share_after_curt     vre_curt vre_share_nocurt
1:                   NA                   NA           NA               NA
2:           0.03597483            0.1967758 0.0001540759        0.1968061
3:           0.01412673            0.1875711 0.0001002543        0.1875899
4:           0.01441916            0.1757875 0.0023678302        0.1762048
5:           0.03032058            0.1547165 0.0453677929        0.1620693
6:           0.06572096            0.1148444 0.1998841529        0.1435347
   total_emissions_co2_milTonnes grid_emissions_factor_kg_kWh
1:                      2065.178                    0.8517769
2:                      1633.143                    0.6735853
3:                      1652.669                    0.6816387
4:                      1676.757                    0.6915735
5:                      1719.900                    0.7093680
6:                      1808.244                    0.7458049
   total_emissions_reduced_milTonnes emissions_reduction
1:                                NA                  NA
2:                          432.0350           0.2091998
3:                          412.5091           0.1997450
4:                          388.4218           0.1880815
5:                          345.2781           0.1671904
6:                          256.9348           0.1244129
   cost_emissions_reduction_pTonneCO2_avg avoided_conv_cap_mw_perVREmw
1:                                     NA                           NA
2:                               7.557954                   0.04509887
3:                               3.108369                   0.10368704
4:                               3.369465                   0.08835494
5:                               7.970634                   0.06808672
6:                              23.216959                   0.04499910
   existing_capacity_coal_MW existing_capacity_ccgt_MW
1:                    189749                     21944
2:                    189749                     21944
3:                    189749                     21944
4:                    189749                     21944
5:                    189749                     21944
6:                    189749                     21944
   existing_capacity_ct_MW  plf_coal    plf_gas
1:                    1579 0.7455022 0.12157051
2:                    1579 0.6359835 0.05515001
3:                    1579 0.6597517 0.06417968
4:                    1579 0.6671672 0.07126608
5:                    1579 0.6672987 0.09520873
6:                    1579 0.6642227 0.13572401

write results csv…

Scenarios for individual plots…

Pull scenario sensitivities…

Plot themes, positions, and labels. Processing 200-400-600 GW scenarios…

Create individual plots [Disabled with eval = FALSE]…

Print and save the plots [Disabled with eval=FALSE]…

Plot functions:
1) Line plots with facets by VRE target
2) Line plots with facets by VRE split

RE share, RE curtailment, Coal PLFs, and grid emissions factors

Wind, solar and coal cost inividual plots

Wind solar and coal Plots in a Panel

Key Insights

Uncurtailed VRE shares range from 17-20% for 200 GW, 30-38% for 400 GW, and 41-50% for 600 GW VRE targets.

However, after curtailment due to technical and economic reasons, VRE shares drop to 11-20% for 200 GW, 12-35% for 400 GW, and 13-43% for 600 GW VRE targets.

VRE curtailment is greatest for solar dominated scenarios. Because of the limit to minimum generation levels (70% in the base scenario), coal plants are unable to back down in order to accommodate solar generation in spit of its zero marginal cost. In this study, coal plants are assumed to stay online throughout the day if their generation is essential to meet net peak demand at any time of the day.

High cost coal, when costs for air pollution control systems are included, decreases average additional cost of VRE to the system to below zero for 3 out of the 5 VRE split scenarios with a 600 GW target because VRE displaces coal generation with a higher variable cost.

Battery storage plots

Modify the battery plots in this code chunk

Minimum Generation Level plots

Modify the min gen plots in this code chunk

Load modified plots

Modify the load modification plots in this code chunk

Additional System Costs due to VRE (and Storage) Panel

Average Emissions Mitigation Costs due to VRE (and Storage) Panel

Changing VRE value with higher hub heights and solar PV tracking

Capacity buildout for base case

Modified load curves